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1.  INTRODUCTION 


The  shipboard  operation  of  optical  laser  systems  is  degraded 

by  turbulence  in  the  marine  planetary  boundary  layer  (MPBL).  The 

2 

refractive  index  structure  function  parameter,  Cn  ,  is  a  measure 

of  the  turbulent  fluctuations  of  refractive  index  over  small 

distances  and  is  the  turbulence  parameter  related  to  optical 

2 

system  performance.  Cn  is  defined  by: 

Cn2  =  [n(x+r )  -  n(x)]2  r“2/3,  (1) 

where  n(x+r)  and  n(x)  are  the  refractive  indices  at  points 
separated  by  a  short  distance  r.  The  overbar  represents  an 
average  value  over  a  spatial  region  or  a  time  interval. 

Refractive  index  fluctuations  cause  a  loss  of  coherence  in  a  laser 
beam,  a  broadening  of  the  beam,  and  a  reduction  in  power  focused 
on  a  target. 

This  climatology  presents  the  means  and  standard  deviations 

O 

of  Cn  values  for  the  North  Atlantic  Ocean  region.  Locales 
characterized  by  large  Cn  values  are  identified  and  indicate 
areas  where  optical  system  performance  will  be  degraded.  The 
diurnal  variation  of  Cn  also  is  presented. 


2.  METHOD 


Surface  layer  turbulence  is  generated  by  the  fluxes  of  heat, 
momentum  and  moisture  in  the  sea  surface  boundary  layer.  The 
values  of  these  fluxes  have  been  related  to  commonly  observed 
meteorological  variables  by  the  development  of  surface  layer 
scaling  laws.  A  model  using  these  laws,  together  with  modifica¬ 
tions  for  use  of  the  model  in  the  MPBL*  has  been  described  by  Burk 
et  al .  (1979).  A  modified  form  of  the  model  is  included  as 
Appendix  A. 

The  meteorological  variables  required  for  the  model  are: 

U  -  wind  speed,  usually  taken  at  10  m  height 

T  -  air  temperature 

T$  -  sea  surface  temperature 

q  -  specific  humidity,  taken  at  10  m  height 
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The  data  base  for  this  study  was  a  subset  of  the  Fleet 
Numerical  Oceanography  Center's  Consolidated  Data  Set  of  marine 
surface  observations.  Eleven  years  of  data  (1963-1973)  were 
extracted  and  processed  for  this  analysis.  This  period  was 
selected  to  avoid  the  inclusion  of  low  precision  merchant  vessel 
data  collected  before  1962.  Husby  (1980)  cites  such  data  as  a 
contributor  to  erroneous  heat  flux  computations. 

The  accuracy  of  the  data  base  can  be  evaluated  by  comparing 
simultaneous  observations  from  merchant  vessels  and  ocean  station 
weather  ships.  Since  merchant  reports  dominate  the  data  base  it 
is  important  that  these  observations  be  reliable.  Bunker  (1976) 
found  good  agreement  between  means  of  wind  speed,  sea  surface 
temperature  and  air  temperature  as  reported  by  merchants  and 
weather  ships  when  at  least  500  observations  are  included  per 
average.  The  uncertainty  in  the  determination  of  the  heat  flux 
can  be  expected  to  be  about  10%  in  such  cases.  This  uncertainty 
increases  as  the  number  of  observations  falls  below  500.  Mean 
values  of  Cn  were  calculated  for  each  1°  latitude  by  1°  longitude 
Marsden  subsquare  in  the  study  area.  Along  the  major  shipping 
lanes  at  least  500  observations  per  month  and  subsquare  generally 
were  reported.  Outside  these  routes  the  number  of  reports 
declined  greatly. 

Several  quality  control  checks  have  been  made  to  improve  the 


reliability  of  the  observations.  Reports  indicating  supersaturated 
conditions  were  discarded  along  with  air  temperatures  outside  the 


range  of  +50°C.  Reported  dew  point  temperatures  were  compared  to 
values  computed  from  reported  wet  and  dry  bulb  temperatures  and 
the  report  discarded  if  the  disagreement  was  >JI°C.  Reported  air- 
sea  temperature  differences  were  similarly  screened.  Sea  surface 
temperatures  <  -8.0°C  were  eliminated  as  were  sea  surface 
temperatures  that  deviated  from  a  running  mean  for  each  month  and 
subsquare  by  more  than  9°C. 

All  observations  were  assumed  to  have  been  made  at  a  height 
of  10  m  since  actual  levels  are  not  reported.  The  effect  of 
differing  heights  is  probably  less  significant  than  other  observa¬ 
tional  errors.  A  minimum  wind  speed  of  1  meter  per  second  was 

allowed.  Reports  of  calm  conditions  were  modified  accordingly. 

2 

Finally,  Cn  was  defined  as  zero  when  the  air  and  sea  surface 
temperatures  were  equal  . 

The  data  successfully  screened  were  grouped  and  analyzed  in 
three  categories,  daytime,  nighttime  and  all  times.  Daytime  was 
defined  to  be  from  0601  to  1800  local  standard  time. 


4 


4.  RESULTS 


2 

The  isopleths  of  mean  Cn  values  are  presented  for  the 
twelve  months  for  day  and  night  (Figure  1-12),  daytime  only 
(Fig.  13-24)  and  for  nighttime  only  (Fig.  25-36).  The  units  are 

_  p  /  O  T  C 

m"  '  x  10  .  Standard  deviations  are  shown  in  Figures  37-72. 

The  most  striking  feature  of  the  figures  is  the  strong 
influence  of  the  Gulf  Stream  on  the  turbulence  values.  The  heat 
flux  from  the  warm  sea  surface  to  the  colder  air  is  accompanied 
by  higher  levels  of  turbulence  than  in  surrounding  areas  of  colder 
water.  The  maximum  values  of  turbulent  intensity  over  the  Gulf 
Stream  occurs  during  winter.  This  corresponds  to  the  strong  heat 
flux  from  the  warm  sea  surface  to  the  cold  winter  air  of 
continental  origin. 

The  diurnal  variation  of  turbulence  appears  to  be  greatest 

during  the  summer  in  the  area  southeast  of  the  Azores.  This  is 

an  area  where  there  is  a  small  annual  variation  of  sea  surface 

temperature  and  wind  speed.  However,  insolation  does  differ 

between  night  and  day,  and  this  contrast  is  largest  during  the 

summer  months.  This  point  leads  us  to  conjecture  that  the 

apparent  diurnal  variation  of  turbulence  levels  between  the  Azores 

and  the  African  coast  is  primarily  due  to  ship  insolation,  and 

2 

does  not  indicate  the  actual  variability  of  Cn  in  the  MPBL.  If 

this  is  correct,  then  the  isopleths  based  only  on  nighttime  data 

most  accurately  represent  the  distribution  of  mean  Cp  values  for 

the  region.  _ 

FIGS.  1-72  ARE  PRESENTED  AS 
A  GROUP  BEGINNING  ON  P.  22. 


APPENDIX  A 


2 

Cn  Turbulence  Model 


This  appendix  describes  a  sequence  of  calculations  to  compute 

2 

the  refractive  index  structure  function  parameter,  Cn  ,  from 
shipboard  measurements  of  air  temperature,  sea  surface  temperature, 
wind  speed  and  specific  humidity.  The  algorithm  is  a  modification 
of  a  preliminary  model  described  by  Burk  et  al .  (1979)  who 

2 

concluded  that  the  surface  layer  scaling  laws  used  to  predict  Cn 
are  accurate  to  a  factor  of  two  and  that  the  accuracy  of  the 
shipboard  observations  introduces  the  same  order  of  uncertainty 

9 

into  Cn  results. 

2 

A . 1  Refractive  index  structure  function  parameter,  Cp 

Since  refractive  index  is  related  to  temperature,  moisture 

and  pressure,  the  fluctuations  of  the  refractive  index  that 

2 

compose  Cn  are  caused  by  fluctuations  in  these  variables.  The 

effect  of  pressure  variations  can  be  neglected  while  variations 

of  temperature  and  moisture  are  described  by  the  temperature 

2 

structure  function  parameter,  CT  ,  the  moisture  structure 

2 

function  parameter,  Cg  ,  and  the  temperature/moisture  correlation 
parameter,  Cgj.  These  component  structure  parameters  are  related 
to  Cn2  by: 


Zn  =  B1CT2  +  B2CeT  +  B3Ce2T 


( A  - 1  ) 


where  B1 ,  Bg  and  are  functions  of  temperature,  pressure  and 
vapor  pressure. 
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Over  the  range  of  normally  observed  values  of  these  structure 

2  2 
parameters,  CT  is  the  largest  contributor  to  C  .  Over  land, 
i  n 

where  moisture  sources  are  relatively  unimportant,  the  contribu- 

2 

tions  of  Ce  and  CeT  are  negligible.  Over  water  the  contributions 

2 

of  these  structure  parameters  are  significant.  To  calculate  Cn 

ir  the  marine  planetary  boundary  layer  (MPBL),  it  is  necessary  to 
2 

calculate  Cy  and  modify  it  for  the  contribution  made  by  moisture 
f 1 uctuati ons . 

In  developing  the  modification  used  in  this  model,  Wesely 

(1976)  expands  the  definition  of  the  refractive  index  structure 

2 

parameter  to  provide  Cp  as  a  function  of  the  component  structure 
functions  and  covariances. 


1  +2r 


eT 


pl-Ag/A^CgTj'  |‘(1-Ag/A1)CeTj 


( A- 2 ) 


where  T  is  the  mean  layer  temperature  in  K,  P  is  the  mean  layer 
pressure  in  mb. 


A]  =  78 . 7 • 1 0~6K  mb 

For  this  analysis  it 
correlated  and  rgy  = 


Ag  =  66.7*10  mb  ^  and  = 

was  assumed  that  e  and  T  are  perfectly 
1 .  Thus , 


C 


2 

n 


( A- 3 ) 


where  B  is  the  Bowen 


ratio  defined  Jjelow. 
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WV4. 


2 

A. 2  Temperature  structure  function  parameter,  Cy 

In  the  surface  layer,  dimensiotal  arguments  first  introduced 

2 

by  Monin  and  Obukhov  can  be  used  to  relate  Cy  to  a  temperature 
scale,  T*;  a  height,  z;  and  a  stability  function,  g(z/L) . 

Wyngaard  et  al .  (1971)  derived  a  semi empi  ri cal  relationship  of 
these  parameters, 

Cy2  =  T*2  z-2/3  g(z/L)  ( A-4 ) 


where  L  is  the  Moni n-Obukho v  length  scale. 

2 

A. 3  Stability  Function  in  Cy 


Using  a  set  of  careful  observations  over  uniform  terrain  in 

Kansas,  Wyngaard  et  al .  (1971)  have  empirically  determined  the 

2 

function  g(z/L)  from  measurements  of  Cy  and  the  scaling  parameters 
L  and  T*.  The  results  have  been  tested  over  water  by  Friehe 
(1977)  and  Davidson  et  al .  (1978).  The  following  expressions  are 
used  in  this  analysis, 


g(z/L)  =  4.9(1 -7z/L)"2/3  (L<0) 

g(z/L)  =  4.9(1+2.75(z/L)2/3)  (L>0). 


( A  -  5  ) 


A . 4  Scaling  functions  in  the  surface  layer 

The  Monin-Obukhov  theory  describes  the  surface  layer  using 
non-dimensional  atmospheric  variables.  Physical  variables  of 
length,  z;  wind  speed,  U;  temperature,  T;  and  specific  humidity,  q; 
are  made  non-dimensional  by  defining  a  scaling  length,  L;  a  scaling 
speed,  U*;  a  scaling  temperature,  T*;  and  a  scaling  humidity,  Q*. 

The  non-dimensional  parameters  describe  the  fluxes  of  momentum, 
heat,  and  humidity  in  the  surface  layer. 
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A . 5  Mon i n-Obukho v  length,  L 

In  a  dry  atmosphere  the  characteristic  length  L^y  is  taken 
as  the  Mon i n -Obukhov  length, 

T  U*2 

L0RY  =  g  kT*  (A‘6) 

where  T  is  mean  surface  layer  temperature  in  K,  g  is  acceleration 
of  gravity,  and  k  is  the  von  Karman  constant.  In  unstable  air 
buoyancy  generated  turbulence  domi nates  at  levels  above  -L,  while 
mechanical  turbulence  dominates  much  below  -L.  The  length  is 
corrected  for  moisture  by  using  the  Bowen  ratio  B, 

L  =  LDRY^  +  1“^  (A-7) 


The  Bowen  ratio  is  related  to  the  heat  and  moisture  fluxes  and  is 
expressed  in  terms  of  the  scaling  temperature  and  humidity. 


8  = 


CpT« 

LWQ* 


A  -  8 ) 


where  C  is  specific  heat  of  air  at  constant  pressure  and  L  is 
P  w 

latent  heat  of  vaporization. 

A. 6  Scaling  wind  speed,  U*;  temperature,  T*;  and  humidity,  Q* 


The  surface  layer  is  characterized  by  fluxes  of  momentum, 
heat  and  moisture  which  are  independent  of  height.  These  constant 
fluxes  are  used  to  define  the  scaling  velocity,  U*;  temperature, 

T+ ;  and  humi di ty,  Q*; 


_  T*  = 

u* 

Q*  = 


T 


(A-9) 


9 


r 


where  U ' ,  T'  and  q'  are  fluctuations  of  horizontal  wind,  tempera¬ 
ture  and  humidity  about  their  means,  and  the  overbar  is  an  average 
among  all  fluctuations.  These  scaling  variables  are  related  to 
the  vertical  profiles  of  each  variable, 

U*  =  kU[ln(^-)  -  ^(z/L)]'1 
o 

Q*  •  |(q-qs)C1n(rj)  -  *q(z/L )]_1  (A-10) 

T*  *  ^•{T-Ts)[ln(7^-)  -  ^(z/L)]"1 

where  U  is  the  wind  speed  in  m/s,  T  and  T$  are  the  air  and  sea 
surface  temperatures  in  °C,  q  and  qg  are  the  specific  humidities 
at  the  observation  height  Z  and  at  the  sea  surface  in  grams  per 
gram.  ZQ,  ZqT,  and  ZQg  are  roughness  parameters  describing 
molecular  transport  of  wind,  temperature,  and  moisture  near  the 
surface,  and  ^M,  ^ ,  and  ifiq  are  empirical  profile  functions.  The 
constant  R  is  included  in  the  equations  for  temperature  and 
humidity.  The  value  of  this  constant  is  accepted  as  0.74  based  on 
reliable  field  experiments. 

A . 7  Profile  functions 

The  profile  functions,  were  formulated  by  Businger  et  al  . 
(1971).  These  functions  have  been  verified  both  over  land  and 
over  water  and  are  given  by. 


10 


1 


*M 


L<0 

■  ln[(l|ii)(l±i)2]-  2  arc  t.„(x)  ♦  ,/t 


*H 


x  "  ^-yn  Z/L)1/4 

*q  ■  2  1"[-^'YHr  ] 


L>0 


♦h  ‘  *6r 

(A- 1 1 ) 


V*Q 


«.  £  L 
R  L 


Constants  in  these  formulations  are  determined  by  fitting  the 
formulations  to  the  data,  and  correspond  best  to  one  particular 
data  set.  Universal  constants  have  so  far  not  been  found. 

The  hypothesis  that  the  ^  are  the  same  for  heat  and  moisture 
is  based  mostly  on  measurements  over  land.  There  is  considerable 
controversy  whether  the  assumption  is  valid  over  the  ocean;  but  no 
better  hypotheses  are  available.  Fortunately,  in  all  equations, 
the  tp 1  s  are  small  compared  to  the  logarithmic  terms,  so  that 
extreme  accuracy  in  the  ^'s  is  not  required.  Further,  ^  occurs 
only  in  the  correction  of  L  for  moisture,  which  itself  is  usually 
a  small  correction. 

A. 8  Roughness  Lengths  ZQ,  ZpT,  Zqq 

The  molecular  transport  of  momentum,  heat  and  moisture  near 
the  surface  is  described  by  the  values  of  ZQ,  ZqT,  and  ZQg 
respectively.  In  the  marine  environment,  the  roughness  lengths 

2 

have  values  on  the  order  of  1  mm  or  less.  The  dependence  of  Cn 

on  ZQ  is  included  primarily  in  the  function  ln^-,  where  Z  is 

0  2 

the  elevation,  taken  to  be  10  meters.  Therefore  Cn  should  be 

insensitive  to  errors  in  Z„ .  Separate  treatments  of  Z  .  lnT  and 

o  o  oT 

Zqq  are  necessary.  The  technique  developed  by  Liu  (1978)  is  used 


m  { 


11 


here.  Liu  derived  relations  between  U  and  U*  and  the  transfer 
coefficients  for  momentum,  heat,  and  moisture,  (Cg,  CH,  Cq). 

These  relationships  are  shown  in  Figure  A-l .  The  roughness  lengths 
can  be  related  to  Cg ,  C^ ,  Cq,  and  the  profile  functions  by, 

lQ  =  Z-exp(J^-  +  *m(z/L)) 

CD 


Z0T  ■  Z-eXP<!?  +  »H<2ZL)) 


zoQ  '  z-e*p<k  EJ-  +  *q<2/L)) 


(A-l 2) 


The  values  of  the  roughness  lengths  cannot  be  obtained  in  closed 

form  expressions,  and  the  following  procedure  was  used: 

30 

a.  An  arbitrary,  initial  value  for  L  of  10  was  chosen. 
This  implies  that  Z/L=0,  (since  Z  was  assumed  to  be  10  m)  and 
indicates  neutral  conditions. 

b.  The  profile  functions,  ij>,  were  assigned  their  neutral 
values  of  zero. 

c.  The  transfer  coefficients  Cq  and  CH  were  given  initial 

_  3 

values  of  1.3  x  10  which  is  approximately  the  maximum  value 
shown  in  Fi gure  A-l . 

d.  The  curve  for  Cq  shown  in  Figure  A-l  was  approximated 
by  the  equation  , 

Cq  =  10‘3-(0. 62375+0. 8414-LogT0(U)) 


(A-l 3) 


Fig.  A-l .  Bulk  transfer  coefficients  at  neutral  conditions  versus 
10  m  wind  U  and  roughness  Reynolds  number  R  (Liu,  1978). 

The  transfer  coefficients  are  approximatelyrcal cul ated  from 
wind  U;  when  scaling  speed  U*  and  roughness  length  Z  are 
found,  transfer  coefficients  are  found  from  R  usi ng°formul as 
in  the  text.  r 


where  U  is  the  wind  speed.  For  this  calculation  only,  reported 
values  of  U < 2 . 5  m/s  were  changed  to  a  minimum  value  of  2.5  m/s  to 
correspond  to  the  lower  limit  of  the  diagram.  Using  Eq.  (A-13), 
an  initial  value  of  Cp  was  computed. 

e.  ZQ,  ZQj  and  ZQq  were  calculated  from  equation  ( A - 1 2 )  . 

A  maximum  value  of  1  mm  was  allowed  for  all  roughness  lengths. 
Larger  values  were  reduced  to  this  limit. 

f.  The  initial  roughness  lengths  were  used  to  calculate  new 
values  for  the  transfer  coefficients,  Cp,  and  Cq.  These  values 
in  turn  were  used  to  compute  the  scaling  parameters  U*,  T*  and 

Q*.  At  this  point  the  roughness  Reynolds  number  was  found  by. 


where  v  is  the  kinematic  viscosity  of  air.  As  shown  in  Figure  A-l 
a  minimum  value  of  0.1  was  allowed  for  Rr. 

g.  A  new  set  of  roughness  lengths  was  calculated  using 
another  approximation  to  the  curve  for  Cp  in  Figure  A-l  and  the 
arbitrary  values  of  =  .0011  and  Cq  =  .0012.  The  equation  for 

CD  is 

CD  =  1 0"3 • ( 1 .233+0. 336-Log10(Rr)).  (A-15) 

_  o 

A  minimum  value  of  0.75x10  was  allowed.  This  corresponds  to  the 
lower  limit  in  Figure  A-l. 


1 

A. 9  Sea  Surface  Temperature,  Ts 

In  most  cases,  bucket  temperature,  T  ,  rather  than  the  surface 

w 

temperature  is  measured.  Tg  was  estimated  from  a  theory  by  Liu 
(1978).  Liu,  (Figure  A-2)  gives  T  -T  as  a  function  of  the  wind 

b  n 

speed  and  the  difference  between  bucket  temperature  and  air 
temperature  at  10  m  height,  aT. 

Liu's  diagram  was  approximated  by  the  following  equations: 

VTs  =  •  025aT  +  0.1  ( U > 6  m  sec"1)  (A-16) 

T  -T  =  O.IaT  +  .8  -  U ( .117  +  . 01 25aT ) ( U<6  m  sec"1) 
w  s 

where  the  temperature  is  in  °C  and  the  wind  speed  is  in  m  sec”1. 

A .10  Evaluation  of  constants 

Certain  empirical  constants  were  introduced  in  the  definition 

9 

of  scaled  quantities,  the  vertical  flux  profiles,  and  in  the  CT 
stability  function.  These  constants  depend  on  the  particular  data 
set  from  which  they  were  determined,  and  there  is  no  uniform 
agreement  in  the  literature  on  the  best  values  of  each.  In 
choosing  the  constants,  a  set  of  values  corresponding  to  one 
particular  data  set  should  be  chosen  because  the  constants  are 
interdependent.  Table  1  provides  the  values  of  the  constants  used 
in  this  analysis. 
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Fig.  A-2.  Temperature  difference  across  the  aqueous  sublayer  versus 
10  m  wind  for  values  of  air-bulk  water  temperature  difference 
Tw"Tio  (Liu,  1978).  These  formulas  account  for  radiation  and 
molecular  effects  in  the  surface  layer. 


Table  1 


Constant 

For  Scaling  Functions 

k 
R 

For  Profile  Functions 

Ym 

yH 
8 

For  Stability  Function 
b 

2 

A. 11  Cal culation  of  C 

_ _ n _ 

?  . 

The  step-by-step  procedure  used  to  calculate  Cn  is  given 
below. 

a.  Modify  the  sea  surface  temperature  measurement  to  an 
actual  sea  surface  temperature  value  using  Eq .  ( A- 1 6 )  and  T  ,  T 

5  W 

and  U . 

30 

b.  Set  initial  stability  to  neutral,  i.e.,  L  =  1 0  ,  Z/iz0, 

and  ip M  =  i^h  =  *Pq  =  °- 

c.  The  bulk  wind  speed  determines  an  initial  value  of  Cq 
using  Eq.  ( A - 1 3 ) .  Initial  values  for  CH  and  Cq  are  set  to  .0013. 

d.  Calculate  ZQ,  ZoT,  and  ZQg  from  Cq,  C^,  Cq  with 
Eq.  ( A- 1 2 ) . 


Value  Used 


0.35 

.74 

15 

6 

4.7 


.667 


e.  Calculate  new  values  for  Cq,  Ch,  and  Cq  by  rearranging 
Eq.  (A-12). 

f.  Evaluate  scaling  parameters  U*,  T*  and  Q*  using  Eq . 

( A  - 1  0 )  . 
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g.  Evaluate  Bowen  ratio,  B,  using  Eq .  (A-8). 

h.  Evaluate  Moni n-Obukhov  length  for  moist  atmosphere  using 
Eqs .  (A-6 )  and  ( A - 7 )  . 

i.  Test  if  the  new  L  is  within  0.1%  of  the  last  value  of 
L.  If  so,  go  to  statement  m.  If  not,  go  to  statement  j. 

j.  Calculate  stability  functions  ^M,  <j>H,  ^  using  Eq.  ( A  - 1 1  ) . 

k.  Calculate  roughness  Reynolds  number  Rr  from  Eq.  ( A- 1 4 ) . 

l.  Evaluate  CQ  from  Eq.  ( A  - 1  5 )  using  Rp  and  set  CH  and  Cg 

to  .0011  and  .0012  respectively.  Go  to  statement  d. 

2 

m.  Evaluate  Cy  stability  function,  g(z/L),  using  Eq .  ( A- 5 ) . 

2 

n.  Evaluate  C-y  using  Eq .  (A-4). 

2 

o.  Calculate  Cn  using  Eq .  (A-3). 

This  procedure  is  graphically  outlined  in  Figure  A-3. 


Figure  A-3.  Flow  chart  of  surface  layer  C  calculation.  Symbols, 
equation  numbers,  and  figure  numbers  arenthose  used  in  text. 


2 

Glossary  and  Values  for  Cp  Turbulence  Model 
Bowen  ratio 

bulk  transfer  coefficients 

moisture  structure  function  parameter 

temperature/moisture  correlation  parameter 

index  of  refraction  structure  function  parameter 

specific  heat  of  air  at  constant  pressure 
(0.24  cal/g  K) 

temperature  structure  function  parameter 

2 

acceleration  of  gravity  (9.8  m/s  ) 
von  Karman  constant  (.35) 

Monin-Obukov  length  (m) 

latent  heat  of  vaporization  (591.7  cal/g) 

optical  index  of  refraction 

atmospheric  pressure  (mb) 

specific  humidity  (g  per  g) 

scaling  humidity  (g  per  g) 

scaling  constant  (.74) 

roughness  Reynolds  number 

temperature  (K  or  °C) 

sea  surface  temperature  (°C) 

bucket  temperature  (°C) 

scaling  temperature  (K) 

wind  speed  (m/sec) 

scaling  velocity  (m/sec) 

observation  height  (assumed  10  m) 

roughness  lengths  (m) 

profile  function  constant  (4.7) 

heat  flux  profile  parameter  (6.0) 

momentum  flux  profile  parameter  (15.0) 

kinematic  viscosity  of  air  (1.35x10”  m  /s) 

profile  functions 
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FIGURES 
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Figures  showing  isopleths  of  mean  values  and  of  standard 
deviations  are  sequenced  numerically  for  the  twelve  months 
of  the  calendar  year  as  follows: 

Means 

Day  and  Night  --  Figs.  1-12  ("All  Times") 

Daytime  Only  --  Figs.  13-24 
Nighttime  Only  -  Figs.  25-36 

Standard  Deviations 

Day  and  Night  --  Figs.  37-48  ("All  Times") 

Daytime  Only  --  Figs.  49-60 
Nighttime  Only  -  Figs.  61-72 
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